Severus: accurate detection and characterization of somatic structural variation in tumor genomes using long reads.

Most current studies rely on short-read sequencing to detect somatic structural variation (SV) in cancer genomes. Long-read sequencing offers the advantage of better mappability and long-range phasing, which results in substantial improvements in germline SV detection. However, current long-read SV detection methods do not generalize well to the analysis of somatic SVs in tumor genomes with complex rearrangements, heterogeneity, and aneuploidy. Here, we present Severus: a method for the accurate detection of different types of somatic SVs using a phased breakpoint graph approach. To benchmark various short- and long-read SV detection methods, we sequenced five tumor/normal cell line pairs with Illumina, Nanopore, and PacBio sequencing platforms; on this benchmark Severus showed the highest F1 scores (harmonic mean of the precision and recall) as compared to long-read and short-read methods. We then applied Severus to three clinical cases of pediatric cancer, demonstrating concordance with known genetic findings as well as revealing clinically relevant cryptic rearrangements missed by standard genomic panels.

A syndromic neurodevelopmental disorder caused by rare variants in PPFIA3.

PPFIA3 encodes the protein-tyrosine phosphatase, receptor-type, F-polypeptide-interacting-protein-alpha-3 (PPFIA3), which is a member of the LAR-protein-tyrosine phosphatase-interacting-protein (liprin) family involved in synapse formation and function, synaptic vesicle transport, and presynaptic active zone assembly. The protein structure and function are evolutionarily well conserved, but human diseases related to PPFIA3 dysfunction are not yet reported in OMIM. Here, we report 20 individuals with rare PPFIA3 variants (19 heterozygous and 1 compound heterozygous) presenting with developmental delay, intellectual disability, hypotonia, dysmorphisms, microcephaly or macrocephaly, autistic features, and epilepsy with reduced penetrance. Seventeen unique PPFIA3 variants were detected in 18 families. To determine the pathogenicity of PPFIA3 variants in vivo, we generated transgenic fruit flies producing either human wild-type (WT) PPFIA3 or five missense variants using GAL4-UAS targeted gene expression systems. In the fly overexpression assays, we found that the PPFIA3 variants in the region encoding the N-terminal coiled-coil domain exhibited stronger phenotypes compared to those affecting the C-terminal region. In the loss-of-function fly assay, we show that the homozygous loss of fly Liprin-α leads to embryonic lethality. This lethality is partially rescued by the expression of human PPFIA3 WT, suggesting human PPFIA3 function is partially conserved in the fly. However, two of the tested variants failed to rescue the lethality at the larval stage and one variant failed to rescue lethality at the adult stage. Altogether, the human and fruit fly data reveal that the rare PPFIA3 variants are dominant-negative loss-of-function alleles that perturb multiple developmental processes and synapse formation.

Monoallelic variation in DHX9, the gene encoding the DExH-box helicase DHX9, underlies neurodevelopment disorders and Charcot-Marie-Tooth disease.

DExD/H-box RNA helicases (DDX/DHX) are encoded by a large paralogous gene family; in a subset of these human helicase genes, pathogenic variation causes neurodevelopmental disorder (NDD) traits and cancer. DHX9 encodes a BRCA1-interacting nuclear helicase regulating transcription, R-loops, and homologous recombination and exhibits the highest mutational constraint of all DDX/DHX paralogs but remains unassociated with disease traits in OMIM. Using exome sequencing and family-based rare-variant analyses, we identified 20 individuals with de novo, ultra-rare, heterozygous missense or loss-of-function (LoF) DHX9 variant alleles. Phenotypes ranged from NDDs to the distal symmetric polyneuropathy axonal Charcot-Marie-Tooth disease (CMT2). Quantitative Human Phenotype Ontology (HPO) analysis demonstrated genotype-phenotype correlations with LoF variants causing mild NDD phenotypes and nuclear localization signal (NLS) missense variants causing severe NDD. We investigated DHX9 variant-associated cellular phenotypes in human cell lines. Whereas wild-type DHX9 was restricted to the nucleus, NLS missense variants abnormally accumulated in the cytoplasm. Fibroblasts from an individual with an NLS variant also showed abnormal cytoplasmic DHX9 accumulation. CMT2-associated missense variants caused aberrant nucleolar DHX9 accumulation, a phenomenon previously associated with cellular stress. Two NDD-associated variants, p.Gly411Glu and p.Arg761Gln, altered DHX9 ATPase activity. The severe NDD-associated variant p.Arg141Gln did not affect DHX9 localization but instead increased R-loop levels and double-stranded DNA breaks. Dhx9-/- mice exhibited hypoactivity in novel environments, tremor, and sensorineural hearing loss. All together, these results establish DHX9 as a critical regulator of mammalian neurodevelopment and neuronal homeostasis.

Hypomyelination caused by a novel homozygous pathogenic variant in FOLR1: complete clinical and radiological recovery with oral folinic acid therapy and review of the literature.

BACKGROUND: Neurodegeneration due to cerebral folate transport deficiency is a rare autosomal recessive disorder caused by biallelic pathogenic variants in FOLR1. Onset typically occurs in late infancy and is characterized by psychomotor regression, epilepsy, and a hypomyelinating leukodystrophy on magnetic resonance imaging. If left untreated, progressive neurodegeneration occurs. However, early treatment with folinic acid has been shown to stabilize or reverse neurological features. Approximately thirty patients have been described worldwide. Here, we report the first two cases with genetically proven cerebral folate transport deficiency from South-Eastern Europe, describe the effect of oral folinic acid therapy on clinical and neuroradiological features and review the literature. RESULTS: Two siblings presented in childhood with clinical and radiological findings consistent with a hypomyelinating leukodystrophy. Exome sequencing revealed a novel homozygous pathogenic variant in FOLR1 (c.465_466delinsTG; p.W156G), confirming the diagnosis of neurodegeneration due to cerebral folate transport deficiency. Folinic acid treatment was promptly initiated in both patients. The younger sibling was treated early in disease course at 2 years of age, and demonstrated complete recovery in clinical and MRI features. The older sibling, who was 8 years of age at the time of diagnosis and treatment, demonstrated partial but substantial improvements. CONCLUSION: We present the first account in the literature that early treatment initiation with oral folinic acid alone can result in complete neurological recovery of both clinical and radiological abnormalities in neurodegeneration due to cerebral folate deficiency. Moreover, through the report of these patients along with review of the literature, we provide information about the natural history of the disease with comparison of treatment effects at different stages of disease progression. This report also reinforces the importance of universal access to genetic testing to ensure prompt diagnoses for treatable disorders.

Rare variants in PPFIA3 cause delayed development, intellectual disability, autism, and epilepsy.

PPFIA3 encodes the Protein-Tyrosine Phosphatase, Receptor-Type, F Polypeptide-Interacting Protein Alpha-3 (PPFIA3), which is a member of the LAR protein-tyrosine phosphatase-interacting protein (liprin) family involved in synaptic vesicle transport and presynaptic active zone assembly. The protein structure and function are well conserved in both invertebrates and vertebrates, but human diseases related to PPFIA3 dysfunction are not yet known. Here, we report 14 individuals with rare mono-allelic PPFIA3 variants presenting with features including developmental delay, intellectual disability, hypotonia, autism, and epilepsy. To determine the pathogenicity of PPFIA3 variants in vivo , we generated transgenic fruit flies expressing either human PPFIA3 wildtype (WT) or variant protein using GAL4-UAS targeted gene expression systems. Ubiquitous expression with Actin-GAL4 showed that the PPFIA3 variants had variable penetrance of pupal lethality, eclosion defects, and anatomical leg defects. Neuronal expression with elav-GAL4 showed that the PPFIA3 variants had seizure-like behaviors, motor defects, and bouton loss at the 3 rd instar larval neuromuscular junction (NMJ). Altogether, in the fly overexpression assays, we found that the PPFIA3 variants in the N-terminal coiled coil domain exhibited stronger phenotypes compared to those in the C-terminal region. In the loss-of-function fly assay, we show that the homozygous loss of fly Liprin- α leads to embryonic lethality. This lethality is partially rescued by the expression of human PPFIA3 WT, suggesting human PPFIA3 protein function is partially conserved in the fly. However, the PPFIA3 variants failed to rescue lethality. Altogether, the human and fruit fly data reveal that the rare PPFIA3 variants are dominant negative loss-of-function alleles that perturb multiple developmental processes and synapse formation.

Cation leak through the ATP1A3 pump causes spasticity and intellectual disability.

ATP1A3 encodes the α3 subunit of the sodium-potassium ATPase, one of two isoforms responsible for powering electrochemical gradients in neurons. Heterozygous pathogenic ATP1A3 variants produce several distinct neurological syndromes, yet the molecular basis for phenotypic variability is unclear. We report a novel recurrent variant, ATP1A3(NM_152296.5):c.2324C>T; p.(Pro775Leu), in nine individuals associated with the primary clinical features of progressive or non-progressive spasticity and developmental delay/intellectual disability. No patients fulfil diagnostic criteria for ATP1A3-associated syndromes, including alternating hemiplegia of childhood, rapid-onset dystonia-parkinsonism or cerebellar ataxia-areflexia-pes cavus-optic atrophy-sensorineural hearing loss (CAPOS), and none were suspected of having an ATP1A3-related disorder. Uniquely among known ATP1A3 variants, P775L causes leakage of sodium ions and protons into the cell, associated with impaired sodium binding/occlusion kinetics favouring states with fewer bound ions. These phenotypic and electrophysiologic studies demonstrate that ATP1A3:c.2324C>T; p.(Pro775Leu) results in mild ATP1A3-related phenotypes resembling complex hereditary spastic paraplegia or idiopathic spastic cerebral palsy. Cation leak provides a molecular explanation for this genotype-phenotype correlation, adding another mechanism to further explain phenotypic variability and highlighting the importance of biophysical properties beyond ion transport rate in ion transport diseases.

IRF7 and UNC93B1 variants in an infant with recurrent herpes simplex virus infection.

Neonatal herpes simplex virus (HSV) infection is a devastating disease with substantial morbidity and mortality. The genetic basis of susceptibility to HSV in neonates remains undefined. We evaluated a male infant with neonatal skin/eye/mouth (SEM) HSV-1 disease, who had complete recovery after acyclovir but developed HSV-1 encephalitis at 1 year of age. An immune workup showed an anergic PBMC cytokine response to TLR3 stimulation but no other TLRs. Exome sequencing identified rare missense variants in IFN-regulatory factor 7 (IRF7) and UNC-93 homolog B1 (UNC93B1). PBMC single-cell RNA-Seq done during childhood revealed decreased expression of several innate immune genes and a repressed TLR3 pathway signature at baseline in several immune cell populations, including CD14 monocytes. Functional studies in fibroblasts and human leukemia monocytic THP1 cells showed that both variants individually suppressed TLR3-driven IRF3 transcriptional activity and the type I IFN response in vitro. Furthermore, fibroblasts expressing the IRF7 and UNC93B1 variants had higher intracellular viral titers with blunting of the type I IFN response upon HSV-1 challenge. This study reports an infant with recurrent HSV-1 disease complicated by encephalitis associated with deleterious variants in the IRF7 and UNC93B1 genes. Our results suggest that TLR3 pathway mutations may predispose neonates to recurrent, severe HSV.

Single-strand mismatch and damage patterns revealed by single-molecule DNA sequencing.

Mutations accumulate in the genome of every cell of the body throughout life, causing cancer and other genetic diseases1-4. Almost all of these mosaic mutations begin as nucleotide mismatches or damage in only one of the two strands of the DNA prior to becoming double-strand mutations if unrepaired or misrepaired5. However, current DNA sequencing technologies cannot resolve these initial single-strand events. Here, we developed a single-molecule, long-read sequencing method that achieves single-molecule fidelity for single-base substitutions when present in either one or both strands of the DNA. It also detects single-strand cytosine deamination events, a common type of DNA damage. We profiled 110 samples from diverse tissues, including from individuals with cancer-predisposition syndromes, and define the first single-strand mismatch and damage signatures. We find correspondences between these single-strand signatures and known double-strand mutational signatures, which resolves the identity of the initiating lesions. Tumors deficient in both mismatch repair and replicative polymerase proofreading show distinct single-strand mismatch patterns compared to samples deficient in only polymerase proofreading. In the mitochondrial genome, our findings support a mutagenic mechanism occurring primarily during replication. Since the double-strand DNA mutations interrogated by prior studies are only the endpoint of the mutation process, our approach to detect the initiating single-strand events at single-molecule resolution will enable new studies of how mutations arise in a variety of contexts, especially in cancer and aging.

X-linked variations in SHROOM4 are implicated in congenital anomalies of the urinary tract and the anorectal, cardiovascular and central nervous systems.

BACKGROUND: SHROOM4 is thought to play an important role in cytoskeletal modification and development of the early nervous system. Previously, single-nucleotide variants (SNVs) or copy number variations (CNVs) in SHROOM4 have been associated with the neurodevelopmental disorder Stocco dos Santos syndrome, but not with congenital anomalies of the urinary tract and the visceral or the cardiovascular system. METHODS: Here, exome sequencing and CNV analyses besides expression studies in zebrafish and mouse and knockdown (KD) experiments using a splice blocking morpholino in zebrafish were performed to study the role of SHROOM4 during embryonic development. RESULTS: In this study, we identified putative disease-causing SNVs and CNVs in SHROOM4 in six individuals from four families with congenital anomalies of the urinary tract and the anorectal, cardiovascular and central nervous systems (CNS). Embryonic mouse and zebrafish expression studies showed Shroom4 expression in the upper and lower urinary tract, the developing cloaca, the heart and the cerebral CNS. KD studies in zebrafish larvae revealed pronephric cysts, anomalies of the cloaca and the heart, decreased eye-to-head ratio and higher mortality compared with controls. These phenotypes could be rescued by co-injection of human wild-type SHROOM4 mRNA and morpholino. CONCLUSION: The identified SNVs and CNVs in affected individuals with congenital anomalies of the urinary tract, the anorectal, the cardiovascular and the central nervous systems, and subsequent embryonic mouse and zebrafish studies suggest SHROOM4 as a developmental gene for different organ systems.

The Role of Dynamic DNA Methylation in Liver Transplant Rejection in Children.

Transcriptional regulation of liver transplant (LT) rejection may reveal novel predictive and therapeutic targets. The purpose of this article is to test the role of differential DNA methylation in children with biopsy-proven acute cellular rejection after LT. Methods: Paired peripheral blood DNA samples were obtained before and after LT from 17 children, including 4 rejectors (Rs) and 13 nonrejectors (NRs), and assayed with MethylC capture sequencing approach covering 5 million CpGs in immune-cell-specific regulatory elements. Differentially methylated CpGs (DMCs) were identified using generalized linear regression models adjusting for sex and age and merged into differentially methylated regions (DMRs) comprising 3 or more DMCs. Results: Contrasting Rs versus NRs, we identified 2238 DMCs in post-LT and 2620 DMCs in pre-LT samples, which clustered in 216 and 282 DMRs, respectively. DMCs associated with R were enriched in enhancers and depleted in promoters. Among DMRs, the proportion of hypomethylated DMRs increased from 61/282 (22%) in pre-LT to 103/216 (48%, P < 0.0001) in post-LT samples. The highest-ranked biological processes enriched in post-LT DMCs were antigen processing and presentation via major histocompatibility complex (MHC) class I, MHC class I complex, and peptide binding (P < 7.92 × 10-17), respectively. Top-ranked DMRs mapped to genes that mediate B-cell receptor signaling (ADAP1) or regulate several immune cells (ARRB2) (P < 3.75 × 10-08). DMRs in MHC class I genes were enriched for single nucleotide polymorphisms (SNPs), which bind transcription factors, affect gene expression and splicing, or alter peptide-binding amino acid sequences. Conclusions: Dynamic methylation in distal regulatory regions reveals known transplant-relevant MHC-dependent rejection pathways and identifies novel loci for future mechanistic evaluations in pediatric transplant subcohorts.

Epimutations in both the TESK2 and MMACHC promoters in the Epi-cblC inherited disorder of intracellular metabolism of vitamin B12.

BACKGROUND: epi-cblC is a recently discovered inherited disorder of intracellular vitamin B12 metabolism associating hematological, neurological, and cardiometabolic outcomes. It is produced by an epimutation at the promoter common to CCDC163P and MMACHC, which results from an aberrant antisense transcription due to splicing mutations in the antisense PRDX1 gene neighboring MMACHC. We studied whether the aberrant transcription produced a second epimutation by encompassing the CpG island of the TESK2 gene neighboring CCDC163P. METHODS: We unraveled the methylome architecture of the CCDC163P-MMACHC CpG island (CpG:33) and the TESK2 CpG island (CpG:51) of 17 epi-cblC cases. We performed an integrative analysis of the DNA methylome profiling, transcriptome reconstruction of RNA-sequencing (RNA-seq), chromatin immunoprecipitation sequencing (ChIP-Seq) of histone H3, and transcription expression of MMACHC and TESK2. RESULTS: The PRDX1 splice mutations and activation of numerous cryptic splice sites produced antisense readthrough transcripts encompassing the bidirectional MMACHC/CCDC163P promoter and the TESK2 promoter, resulting in the silencing of both the MMACHC and TESK2 genes through the deposition of SETD2-dependent H3K36me3 marks and the generation of epimutations in the CpG islands of the two promoters. CONCLUSIONS: The antisense readthrough transcription of the mutated PRDX1 produces an epigenetic silencing of MMACHC and TESK2. We propose using the term 'epi-digenism' to define this epigenetic disorder that affects two genes. Epi-cblC is an entity that differs from cblC. Indeed, the PRDX1 and TESK2 altered expressions are observed in epi-cblC but not in cblC, suggesting further evaluating the potential consequences on cancer risk and spermatogenesis.

Genetic perturbation of PU.1 binding and chromatin looping at neutrophil enhancers associates with autoimmune disease.

Neutrophils play fundamental roles in innate immune response, shape adaptive immunity, and are a potentially causal cell type underpinning genetic associations with immune system traits and diseases. Here, we profile the binding of myeloid master regulator PU.1 in primary neutrophils across nearly a hundred volunteers. We show that variants associated with differential PU.1 binding underlie genetically-driven differences in cell count and susceptibility to autoimmune and inflammatory diseases. We integrate these results with other multi-individual genomic readouts, revealing coordinated effects of PU.1 binding variants on the local chromatin state, enhancer-promoter contacts and downstream gene expression, and providing a functional interpretation for 27 genes underlying immune traits. Collectively, these results demonstrate the functional role of PU.1 and its target enhancers in neutrophil transcriptional control and immune disease susceptibility.

Single-cell analysis of human adipose tissue identifies depot and disease specific cell types.

The complex relationship between metabolic disease risk and body fat distribution in humans involves cellular characteristics which are specific to body fat compartments. Here we show depot-specific differences in the stromal vascual fraction of visceral and subcutaneous adipose tissue by performing single-cell RNA sequencing of tissue specimen from obese individuals. We characterize multiple immune cells, endothelial cells, fibroblasts, adipose and hematopoietic stem cell progenitors. Subpopulations of adipose-resident immune cells are metabolically active and associated with metabolic disease status and those include a population of potential dysfunctional CD8+ T cells expressing metallothioneins. We identify multiple types of adipocyte progenitors that are common across depots, including a subtype enriched in individuals with type 2 diabetes. Depot-specific analysis reveals a class of adipocyte progenitors unique to visceral adipose tissue, which shares common features with beige preadipocytes. Our human single-cell transcriptome atlas across fat depots provides a resource to dissect functional genomics of metabolic disease.

Computational Analysis of HLA-presentation of Non-synonymous Recipient Mismatches Indicates Effect on the Risk of Chronic Graft-vs.-Host Disease After Allogeneic HSCT.

Genetic mismatches in protein coding genes between allogeneic hematopoietic stem cell transplantation (allo-HSCT) recipient and donor can elicit an alloimmunity response via peptides presented by the recipient HLA receptors as minor histocompatibility antigens (mHAs). While the impact of individual mHAs on allo-HSCT outcome such as graft-vs.-host and graft-vs.-leukemia effects has been demonstrated, it is likely that established mHAs constitute only a small fraction of all immunogenic non-synonymous variants. In the present study, we have analyzed the genetic mismatching in 157 exome-sequenced sibling allo-HSCT pairs to evaluate the significance of polymorphic HLA class I associated peptides on clinical outcome. We applied computational mismatch estimation approaches based on experimentally verified HLA ligands available in public repositories, published mHAs, and predicted HLA-peptide affinites, and analyzed their associations with chronic graft-vs.-host disease (cGvHD) grades. We found that higher estimated recipient mismatching consistently increased the risk of severe cGvHD, suggesting that HLA-presented mismatching influences the likelihood of long-term complications in the patient. Furthermore, computational approaches focusing on estimation of HLA-presentation instead of all non-synonymous mismatches indiscriminately may be beneficial for analysis sensitivity and could help identify novel mHAs.

H3K27M induces defective chromatin spread of PRC2-mediated repressive H3K27me2/me3 and is essential for glioma tumorigenesis.

Lys-27-Met mutations in histone 3 genes (H3K27M) characterize a subgroup of deadly gliomas and decrease genome-wide H3K27 trimethylation. Here we use primary H3K27M tumor lines and isogenic CRISPR-edited controls to assess H3K27M effects in vitro and in vivo. We find that whereas H3K27me3 and H3K27me2 are normally deposited by PRC2 across broad regions, their deposition is severely reduced in H3.3K27M cells. H3K27me3 is unable to spread from large unmethylated CpG islands, while H3K27me2 can be deposited outside these PRC2 high-affinity sites but to levels corresponding to H3K27me3 deposition in wild-type cells. Our findings indicate that PRC2 recruitment and propagation on chromatin are seemingly unaffected by K27M, which mostly impairs spread of the repressive marks it catalyzes, especially H3K27me3. Genome-wide loss of H3K27me3 and me2 deposition has limited transcriptomic consequences, preferentially affecting lowly-expressed genes regulating neurogenesis. Removal of H3K27M restores H3K27me2/me3 spread, impairs cell proliferation, and completely abolishes their capacity to form tumors in mice.

Very long intergenic non-coding RNA transcripts and expression profiles are associated to specific childhood acute lymphoblastic leukemia subtypes.

Very long intergenic non-coding RNAs (vlincRNAs) are a novel class of long transcripts (~50 kb to 1 Mb) with cell type- or cancer-specific expression. We report the discovery and characterization of 256 vlincRNAs from a cohort of 64 primary childhood pre-B and pre-T acute lymphoblastic leukemia (cALL) samples, of which 61% are novel and specifically expressed in cALL. Validation was performed in 35 pre-B and pre-T cALL primary samples. We show that their expression is cALL immunophenotype and molecular subtype-specific and correlated with epigenetic modifications on their promoters, much like protein-coding genes. While the biological functions of these vlincRNAs are still unknown, our results suggest they could play a role in cALL etiology or progression.

Promoter capture Hi-C-based identification of recurrent noncoding mutations in colorectal cancer.

Efforts are being directed to systematically analyze the non-coding regions of the genome for cancer-driving mutations1-6. cis-regulatory elements (CREs) represent a highly enriched subset of the non-coding regions of the genome in which to search for such mutations. Here we use high-throughput chromosome conformation capture techniques (Hi-C) for 19,023 promoter fragments to catalog the regulatory landscape of colorectal cancer in cell lines, mapping CREs and integrating these with whole-genome sequence and expression data from The Cancer Genome Atlas7,8. We identify a recurrently mutated CRE interacting with the ETV1 promoter affecting gene expression. ETV1 expression influences cell viability and is associated with patient survival. We further refine our understanding of the regulatory effects of copy-number variations, showing that RASL11A is targeted by a previously identified enhancer amplification1. This study reveals new insights into the complex genetic alterations driving tumor development, providing a paradigm for employing chromosome conformation capture to decipher non-coding CREs relevant to cancer biology.

Genome-wide association study implicates immune dysfunction in the development of Hodgkin lymphoma.

To further our understanding of inherited susceptibility to Hodgkin lymphoma (HL), we performed a meta-analysis of 7 genome-wide association studies totaling 5325 HL cases and 22 423 control patients. We identify 5 new HL risk loci at 6p21.31 (rs649775; P = 2.11 × 10-10), 6q23.3 (rs1002658; P = 2.97 × 10-8), 11q23.1 (rs7111520; P = 1.44 × 10-11), 16p11.2 (rs6565176; P = 4.00 × 10-8), and 20q13.12 (rs2425752; P = 2.01 × 10-8). Integration of gene expression, histone modification, and in situ promoter capture Hi-C data at the 5 new and 13 known risk loci implicates dysfunction of the germinal center reaction, disrupted T-cell differentiation and function, and constitutive NF-κB activation as mechanisms of predisposition. These data provide further insights into the genetic susceptibility and biology of HL.

Hidden genomic MHC disparity between HLA-matched sibling pairs in hematopoietic stem cell transplantation.

Matching classical HLA alleles between donor and recipient is an important factor in avoiding adverse immunological effects in HSCT. Siblings with no differences in HLA alleles, either due to identical-by-state or identical-by-descent status, are considered to be optimal donors. We carried out a retrospective genomic sequence and SNP analysis of 336 fully HLA-A, -B, -DRB1 matched and 14 partially HLA-matched sibling HSCT pairs to determine the level of undetected mismatching within the MHC segment as well as to map their recombination sites. The genomic sequence of 34 genes locating in the MHC region revealed allelic mismatching at 1 to 8 additional genes in partially HLA-matched pairs. Also, fully matched pairs were found to have mismatching either at HLA-DPB1 or at non-HLA region within the MHC segment. Altogether, 3.9% of fully HLA-matched HSCT pairs had large genomic mismatching in the MHC segment. Recombination sites mapped to certain restricted locations. The number of mismatched nucleotides correlated with the risk of GvHD supporting the central role of full HLA matching in HSCT. High-density genome analysis revealed that fully HLA-matched siblings may not have identical MHC segments and even single allelic mismatching at any classical HLA gene often implies larger genomic differences along MHC.